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Recently a model has been proposed that links dark matter and neutrino masses. The dark 
matter candidate which is dubbed as SLIM has a mass of MeV scale and can show up at low energy 
experiments. The model also has a high energy sector which consists of a scalar doublet, (<f>~ , (f>°). 
' We discuss the potential of the LHC for discovering the new scalars. We focus on the 4> + 4>~ and 

y—i ■ 4' ±c t >0 production and the subsequent decay of the charged scalar to a charged lepton and the SLIM 

which appears as missing energy. Identifying the background, we estimate the signal significance and 
£S) ' find that it can exceed 5a at 30 fb" 1 for the 14 TeV run at the LHC. We comment on the possibility 

_ , of extracting the flavor structure of the Yukawa couplings which also determine the neutrino mass 

matrix. Finally, we discuss the prospects of this search at the current 7 TeV run of the LHC. 
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<—i I. INTRODUCTION 

The nature of dark matter and tiny but nonzero values of neutrino masses are two unsolved mysteries in particle 
physics and cosmology in early 21st century. Various models have been developed that explain nonzero neutrino 
masses. Rich literature also exists on the models suggesting a dark model candidate. Attempts to link these two 
i—i, mysteries and to solve them within a unique framework have only recently been started. In particular, in Ref. [l| a 
simple low energy scenario has been introduced which explains both phenomena in a minimalistic way. The scenario 
adds a real scalar which plays the role of the dark matter as well as two (or more) right-handed Majorana neutrinos that 
along with the scalar couple to the left-handed neutrinos. Through this Yukawa coupling, the left-handed neutrinos 
^vq | acquire a small Majorana mass at one-loop level. A Zi symmetry under which all new particles are odd is imposed 
to guarantee the stability of the lightest new particle which plays the role of the dark matter and is dubbed as SLIM. 
The Z2 symmetry forbids a Dirac mass term for neutrinos. Within this scenario the main annihilation mode for dark 
, matter is to (anti-)neutrino pairs. To simultaneously satisfy upper bound on the neutrino masses and to account for 
"5 ' the observed dark matter abundance within thermal scenario (i.e., (c(DM + DM anything)?;) ~ 3 ■ 1CT 26 cm 3 /sec), 
at least one of the right-handed neutrinos as well as the scalar have to be lighter than 10 MeV. On the other hand, 
• these conditions lead to a lower bound on the couplings which makes the scenario testable in low energy experiments 
^! I such as the measurement of Br(7r, K — > l a + missing energy) @. 
■ ■ *h , With this minimal content, the new Yukawa coupling is not invariant under the electroweak symmetry so the low 
energy effective scenario has to be embedded in a SU(3)xSU(2)xU(l) invariant model at higher energies. An example 
of such a model is introduced in Q . The model includes a scalar electroweak doublet which has Yukawa couplings with 
charged leptons and the right-handed neutrinos. Like in the case of scenario [l[, this Yukawa coupling gives rise to the 
neutrino masses at one-loop level. The model inherits the testability at low energy experiments from the low energy 
scenario that has been embedded in it. In addition, observable effects are expected in high energy experiments such 
as the LHC as well as in searches for lepton flavor violating decays such as fi — > e-f. If the masses of the components 
of the electroweak doublet are not too high, these scalar particles can be produced through electroweak interactions 
at the LHC. As indicated in Ref. [3], it is in principle possible to derive the flavor structure of the Yukawa couplings 
at the LHC and compare it with the flavor structure deduced from the neutrino mass matrix. 

In the present paper, we elaborate on this idea in detail. Considering various backgrounds and by employing 
state-of-the-art techniques for reconstructing the signal, the possibility of discovering the new particles is studied. We 
then discuss and propose conditions under which the Yukawa couplings can be derived. Several packages are linked 
together to ensure a reasonable framework for the simulation and the analysis of events. 

The paper is organized as follows. In sect. HU a brief description of the model is presented. In sect. IIII1 the signal 
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and background events are identified and tools for establishing the simulation and analysis framework are introduced. 
The signal and background cross sections are then calculated. In sect. HVl the number of expected events at LHC 
collected with a beam of 14 TeV center of mass energy after 30 integrated luminosity is calculated. Full set of 
kinematic cuts are designed step by step and their efficiencies on both signal and background events are calculated 
and shown in tables. The signal significance for each studied category and each final state is calculated. In sect. 
IVl a brief discussion on the possibility of measuring the model parameters is presented. In sect. I VI) an alternative 
discovery channel is presented. In sect. IVII1 the discovery potential of the 7 TeV run of the LHC is discussed. The 
paper ends with a conclusion section on the discovery potential of the model and its new particles. 



II. THE MODEL 

The model has a minimalistic content and adds only the following fields to the Standard Model: 

• Two (or more) right-handed neutrinos, Ni\ 

• A scalar SU(2) doublet $ which carries hypercharge: $ = (</>° <j>~) where (f>° — ((f>i + i(f>2)/V2 with real <pi and 

• A singlet scalar rj. 

A Z2 symmetry is imposed on the model under which these new particles are odd but the SM particles are even. The 
Z2 symmetry forbids the Yukawa coupling of form • L and makes the lightest particle stable. The general Z2 
invariant renormalizable Lagrangian involving the scalars can be written as 
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We are interested in the regime that only the standard Higgs, acquires a vacuum expectation value so the Z2 
symmetry is maintained. After the electroweak symmetry breaking, the third term in Lagrangian ([T]) leads to a 
mixing between r\ and the neutral component of For simplicity, we take the CP invariant case so will be real 
and as a result, r\ will mix only with the CP-even component of </>°. The mass eigenstates 5% and 62 can be written as 
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The CP-odd component of $ remains mass eigenstate with mass eigenvalue 

2 2 2 

™% =m% + \i V -f + \ 1 V -f-\ 2 V -f. 
The mass of the charged component of $ is 



4- =™4 + A 4^f. 



Terms involving right-handed neutrinos are 



C = -g ia Ni& ■ L a 



(6) 



where L a is the lepton doublet of flavor a: = (i/£„ £ La ). As shown in Q, this coupling at one-loop level leads to 
the neutrino mass matrix as follows 
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With only one right-handed neutrino, the neutrino mass matrix will have two zero mass eigenvalues and cannot 
accommodate the neutrino data. The model has to include at least two right-handed neutrinos. From Eq. ([?)) ■ it 
is straightforward to show that with only two right-handed neutrinos, the determinant of the neutrino mass matrix 
vanishes which means one of the neutrino mass eigenvalues is zero and the neutrino mass scheme is hierarchical. By 
adding another right-handed neutrino non-hierarchical neutrino mass scheme can also be obtained. Here, we however 
concentrate on the most economic case with only two right-handed neutrinos. In the appendix, we derive constraints on 
the flavor structure of the couplings from the neutrino mass matrix. As shown in the appendix, in the case of inverted 
hierarchical scheme |<7j^| ~ \9ir\- Two specific solutions leading to normal hierarchical mass schemes for neutrinos are 
shown in Tab. Q] To obtain these results, we have set Am^ = 8x 10~ 5 eV 2 and Am 2 tm = 2.5 x 10~ 3 eV 2 and 
A2 = which leads to = m^, 2 . N 2 is taken to be heavy and the mixing between 77 and <f>\ is taken to be 0.01. 6 
is the arbitrary mixing angle that appears in the coupling structure (see Eqs. (|AllA2p V Notice that by varying 
and A2 (i.e., parameters determining the heavy sector masses), the flavor structure (i.e., gia/gip) does not change. 
Taking m^ )1 — — 150 GeV, g\ a is derived from Eq. (IA1[) in the appendix and shown in the last line of the table. 

In this model, Si, being the lightest new particle and therefore stable, plays the role of the dark matter. We shall 
call S\ "SLIM" . Dark matter abundance is fixed by 

2 

= 3 • 10~ 26 cm 3 /sec . (8) 

As shown in |3j, by combining the information on a tot from the dark matter abundance with m„, we find 

ms 1 < rriN 1 ~ few MeV . 

That is while the other scalar particles, S2, 4>2 and <p~ have to be heavy enough to avoid the direct bounds from 
accelerators. The strongest model-independent direct bound on the charged Higgs is 80 GeV [4] and the bound 
on ms 2 and m^ 2 is about 90 GeV (5|. At first sight, it might seem counterintuitive that while W~ and (f>~ with 
m^- = 80 GeV have the same mass and charge, the former is already discovered but such a <p~ has escaped detection. 
This can be explained by the difference in their spin: In the process // — > 7* — > (f)~4> + , the angular momentum 
conservation implies that the final particles should be produced in the p-wave which gives rise to a suppression by an 
extra factor of (1 — 4m 2 _/s) which for the maximum LEP energy amounts to 0.4. Moreover, the W~W + production 
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TABLE I: Model parameters. 



can take place via the £-channel neutrino exchange which has no counterpart in the <p + 4>~ production. To our best 
knowledge, no analysis has been carried out on the Tevatron data that applies to our model. In principle, valence 
quarks and antiquarks inside the proton and antiproton beams at Tevatron can give rise to the 4> + 4>^ production 
via the s-channcl Z* ,7* exchanges. However, each of the valence quarks on average carries only about 0.16 of the 
proton energy so the effective center of mass energy will be around 300 GeV(~ 0.16 x 2 TeV). As a result like the 
LEP case, the production rate at Tevatron will be suppressed. The same process also plays the dominant role in 
the <p + <j)~ production at the LHC but in the case of the LHC, the antiquark taking part in the production will be 
a sea quark which on average carries a lower fraction of the proton energy, (xq) ~ 0.027. The average center of 
mass energy of quark antiquark pair for the 14 TeV run will be */ (%q) (%q)E cm ~ 800 GeV which is still higher than 
that for the Tevatron despite the fact that in the case of the Tevatron, the quark antiquark pair are both valence 
quarks but in the case of the LHC one of the quarks is a sea quark. We shall discuss the effect in detail. As we 
shall demonstrate the discovery chance rapidly increases with center of mass energy so the discovery potential of the 
Tevatron cannot be significant comparing to that of the LHC. In addition to quark antiquark annihilation mode, 
4> + <p~ can be produced via gluon fusion. Similar consideration holds for the production of 82 and </>2- As we shall 
see, the suppression discussed above has the negative consequence for our analysis: The signal for the 4>~ production 
suffers from large background from the W~ production. 

This model has similarities with the so-called inert model [6] but here we have an extra singlet scalar and the main 
annihilation mode of dark matter pair is into neutrinos. Nevertheless, the contribution to the oblique parameters in 
our model is similar to that in the inert model. Like the inert model, the new contribution to the oblique parameters 
can cancel the one from a SM Higgs so heavier masses for Higgs can be made compatible with the SM. Throughout 
this paper, we assume that rru- is equal or lighter than mj 2 and m^ 2 . Moreover, we shall assume that m<f, 2 — TOa-, 
ms 2 — m<f,- and \mg 2 — m^ 2 \ do not exceed 80 GeV to forbid two body decays #2,02 — > W + 4 } '~ or <f>2 — > 82Z (or 
8 2 —¥ 4>2Z) ■ This assumption is just for simplicity of the analysis. Probing the whole parameter space is beyond the 
scope of the present paper. 

In this model, 4>~ dominantly decays via the gi a couplings: 

| fl , I 2 (m 2 - mlr ) 2 

r(0" l a N % ) = 7°' 3 NJ for m Ni < m,- . (9) 
167r mr,_ 
4> 

As discussed earlier, N\ has to be light so decay modes <fr — > A^e - , and N\t~ are guaranteed to exist. If the 

mass of the other right-handed neutrino, Nk, is lighter than mj,- but mjv t ~ fn<h-i by studying the energy spectrum 
of the emitted charged lepton, one can in principle distinguish between this decay and (j>~ — > i^Ni. Moreover, M^ k 
and |<7fe Q | can in principle be derived. If for all k larger than one mjv fc > m$- , the decay mode 4>~ — ► ^fc^a ^ s n °t °P en 
so the decay mode — s- <^+missing energy) gives only |gi a | 2 . In case that rriN k ~ WiVu ^ w iU n °t in practice be 
possible to distinguish between N\ and Nk so the signal —> £~+ missing energy) will practically yield ^ \gia\ 2 
where the sum runs over all Ni for which mjv 4 < m*-. In this paper, we focus on the case with only two with 
mjv 2 > JTU-. Apart from this decay mode, 4>~ can have other decay modes such as (f>~ — > W~5\, (\T — > 8\t^v and 

— > 8\W^^ which as seen from Fig. Q] are all subdominant for values of g\ a ~ 0.01 and sin a ~ 0.01. 
The hypercharge of the new doublet, $, is the same as that of Hd in the MSSM. Thus, in the limit that (Hd) 
vanishes (i.e., tan/3 — > 0) and the r]-(f> mixing goes to zero (i.e., a — > 0), the gauge interactions of <f>2 and <jr^ are 
similar respectively to those of A and H ± . Depending on the value of m ( p 1 , 8 2 will have similar gauge interactions 
as any of the CP-even Higgses, h or H within the MSSM. The pair production of these particles at the LHC takes 
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FIG. 1: Branching ratio of <j> decays as a function of its mass. Point A from Tab. [Jhas been used as input for this calculation. 

place through the electroweak interactions. That is the <j) + <j)~ , (fi^fo and (^82 production are respectively similar to 
H + H~ , H^AP and H^h(H) production within the MSSM with tan/3 — > 00. However, the Yukawa couplings of $ in 
our model is quite different from those of within the MSSM. As a result, the decay processes for the $ components 
differ from those for the Higgs components. Neutral CP-odd Higgs within the MSSM can be singly produced via 
gluon fusion but the parallel does not exist within this model in which Yukawa coupling between $ and the b quarks 
are forbidden. 

To perform our analysis, we focus on the points shown in Tab. HI The point A in Table Q] is chosen to search for 
tauonic and muonic final states while the point B is used to search for final states involving electrons. The point A 
in the parameter space leads to a negligible branching ratio for <p~ decay to electrons while decay to a fi or r lepton 
is possible with equal probability. That is at point A, BR(0 ± —> /i ± iVi) = BR(0 ± — > t Ni) ~ 0.5. Figure [T] shows 
the contributions from all the decay channels and verifies that the above approximation is valid up to a maximum 
1% error for m^i) = 130 GeV. 

In this scenario three samples of events are expected corresponding to different final states: fifj,, fir and tt. The 
tt final state suffers from the large hadronic backgrounds such as W+jets with W decaying to r or jets. These 
background events may survive due to light jets faking the r. 

As we shall see, the fi/j, and fir final states can be observed with a significance exceeding 5a at 30/6 -1 . The next 
sections are devoted to signal and background simulation. The event analysis and selection is then described based on 
the kinematic cuts applied in order to increase the signal to background ratio. A full study of kinematic distributions 
is performed and a dedicated t jet identification is applied similar to the standard LHC algorithms in Ref. 0- Finally 
the signal significance is estimated for each search category and final state. 



III. SIGNAL AND BACKGROUND SIMULATION 



Considering the above observation the ± pair production is generated using the process pp — > H + H~ setting 
tan/3 = 100 which is reasonably high. These events are generated using PYTHIA 8.1.2.5 @. To account for the 
t leptons proper decays, TAUOLA package [|| is used. Since PYTHIA 8 is based on C++ programming language, 
TAUOLA C++ interface [13] is used to link TAUOLA and PYTHIA 8 for the correct production of t hadronic decays. 
Since TAUOLA interface output is in the HepMC format, HepMC 2.5.1 [TT[ is linked to TAUOLA interface and is 
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FIG. 2: Signal cross section as a function of m^i . 



Process 


W+W~ 


tt 


W+jets 


Z+jets 


Cross Section 


115.5±0.4 pb 


878.7±0.5 pb 


187.1±0.1 nb 


258.9±0.7 nb 



TABLE II: Background cross sections calculated using MCFM package. 



used in the analysis. 

Throughout the analysis, the r leptons are identified through their hadronic decays which results in a narrow r jet 
which will be described in detail later in the next sections. 

The jet reconstruction is performed using FASTJET 2.4.1 [l2[ with anti-kt algorithm [l3| and a cone size of 0.4 and 
the Et recombination scheme. 

The signal cross section is calculated using PYTHIA. As mentioned earlier, the charged Higgs pair production is 
expected to simulate ± events if a large enough tan/3 is used. Since the signal cross section decreases rapidly with 
increasing m^± , the low mass (fr^ is studied in this analysis. Figure [5] shows the signal cross section as a function of 
m^i. The main background events are W + W~, tt, W+jets and Z+jets. Their cross sections are calculated using 
MCFM 5.7 [lij at NLO. Table [TT1 shows the background cross sections calculated using MCFM with renormalization 
and factorization scale set to the top quark mass (172.5 GeV) for tt events, the W boson mass (80.4 GeV) for W + W~ 
and VF+jets events and the Z boson mass (91 GeV) for Z+jets events. For both event simulation and cross section 
calculation purposes, MRST 2004 NNLO PDF is used by linking LHAPDF 5.8.1 QJ] to the PYTHIA event generator 
as well as MCFM. 

In order to calculate the expected number of events at 30/6 -1 , the total cross sections in Table [TT1 are multiplied by 
the proper branching ratios. For tt events, while one of W bosons (from i or t decay) decays to a muon, the other 
W, is forced to decay to a r lepton or light jets. The same decay channel is considered for W + W~ events. The W 
boson decay to light jets is considered to account for the r jet fake rate from light jets. For VF+jet events, the W 
boson is forced to decay to a muon and appearance of a r jet is only possible again through the r jet fake rate. The 
Z+jets process is only considered for the final state and the Z boson is forced to decay to a pair of muons. In this 
analysis the following values are used for branching ratio of W and Z bosons according to Particle Data Group [l6| : 
BR(W+ fi+v) = 10.57%, BR(T^+ ->■ t+v) = 11.25%, BR(W+ hadrons) = 67.60%, BR(Z ->■ fj.fi) = 3.36%. 
As is seen from Table [TT1 a large suppression factor is needed due to the large background cross sections. In the 
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next sections, it is shown that a reasonable background suppression can be achieved by applying kinematic cuts and 
using a state-of-the-art r identification and selection. 



IV. EVENT KINEMATICS 



As mentioned earlier in this analysis only two final states of and r/z from (j) + 4>~ decay are considered. In the 
following, two separate analyses are designed and described in detail. It is convenient to define the transverse energy, 
Et, as 




where the z direction is taken parallel to the beam direction. 



A. The t fi ET as Final State 



This final state is searched for by requiring one muon and one r jet. The muon transverse momentum distribution 
has been shown in Fig. [3] The small contribution from soft muons in r lepton decays (see Fig. |4|), can be viewed 
as a small excess in signal and background events. In the signal events, in the search for r /i E™' lss final state, this 
contribution may come from pp — > </> + </>~ — > ttE™ iss , with one r decaying hadronically and the other decaying to a 
muon. However, due to the small branching ratio of r decay to muons (~ 17%) and the soft pr distribution of these 
muons, no sizable contribution from these events is expected. The same considerations hold for muonic decays of the 
r leptons in background events like ti and WW. As a result, the r leptons mainly contribute in the analysis through 
their hadronic decays. The following kinematic requirement is applied on muons in the event: 



p£ > 50 GeV, \t}\ < 2.5 



(10) 



where r\ is the pseudorapidity defined as 



rj = — log 



tan - 
2 



in which 9 is the angle between the momentum and the beam axis, 
backgrounds with soft muons in the final state such as W+jets. 



This kinematic cut is applied mainly to suppress 



The r jet identification is performed by first reconstructing jets and requiring kinematic cuts as the following: 

E^ > 30 GeV, \r)\ < 2.5 . (11) 

A reconstructed jet which passes the above requirement is considered as a r jet candidate. All selected jets should be 
separated enough from muons with px > 20 GeV with the following requirement: 

Ai? (joUl) > 0.4 . (12) 

Here AR is defined as AR = y/ (Aij) 2 + (A(p) 2 where r\ is the pseudorapidity and 4> is the azimuthal angle with 
respect to the collider beam axis. Figure [5] shows the transverse energy distribution of selected jets and the number of 
reconstructed jets passing above requirements are shown in Figure [6] Due to the requirements of Eq. [TT]and[T2l some 
signal events are killed, however, a larger fraction of background events fall in zero jet bin. The single W production 
has the largest cross section when produced with no jet so a large fraction of them fill the zero jet bin. Another 
reason is the fact that the jet accompanying the W boson is usually a spectator quark flying in the forward-backward 
direction and as a result, failing the requirement of Eq. 1111 The r hadronic decay produces predominantly one or 
three charged pions with a sizably high transverse momentum due to the low charged track multiplicity compared 
to the light quark jets. The leading track associated to the jet is found by searching for the highest-pr track in the 
matching cone defined around the jet axis with an opening angle of AR < 0.1. 

The desired polarization of the r leptons in ± decays is not simulated by using charged Higgs events in PYTHIA 
because of opposite polarizations between the PYTHIA charged Higgs events and our model. Thus, in this analysis 
polarization dependent cuts on leading track pp, the ratio of leading track pj< over r jet energy, and EJp lSB are removed. 
The TAUOLA package is still needed to get correct branching ratio of r lepton decays and produce hadronic decays 
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FIG. 3: Muon transverse momentum distributions in signal and background events in the r/i_BJ" ss final state. 



with proper charged track multiplicity and decay topology. The rest of r identification algorithm is applied as the 
following. The isolation cone is defined with a cone size of AR < 0.4 around the leading track. The signal cone is 
also defined around the leading track with AR < 0.07. The isolation is applied by requiring no charged track with 
Pt > 1 GeV to be in the isolation annulus defined as 0.07 < AR < 0.4. This requirement suppresses light quark jets 
which contain many charged tracks in their isolation annulus. 

To further suppress the fake r jet rate and increase the r jet purity in events with genuine r jets, a search for charged 
tracks in the signal cone is made and the number of tracks in the signal cone is required to satisfy the following 
requirement: 

Number of signal tracks — 1 or 3 . (13) 

Figure [7] shows distribution of the number of signal tracks in signal events before applying the cut. As seen from Fig. 
[3 r jets have undergone 1- or 3-prong decays predominantly. The muons and r leptons are expected to be produced 
back-to-back in the signal events. Figure |5] shows the distribution of the azimuthal angle between the two leptons. In 
the case of the r leptons, the r jet direction is used. As is observed from Fig. 13 background events also show similar 
distributions. Thus, no cut on A(f> between the muon and the r jet is applied. In the next step, the r lepton charge is 
calculated as the sum of charges of tracks in the signal cone. Since muons and r leptons in signal events are produced 
with opposite charges, the following requirement is applied: 

Muon charge + r jet charge = . (14) 

The analysis is therefore only sensitive to the sum of charges of the two leptons and does not care whether /j, and r 
are individually positive or negative. 

An event is selected if it passes all the above selection requirements. Table Hill shows selection efficiencies for signal 
events with TO(0±) ranging from 80 GeV to 130 GeV. The quoted efficiencies show fraction of events which pass a cut 
and contain exactly the same number of physical objects as the cut requires, e.g., the fourth row shows the fraction 
of events having exactly one muon passing the kinematic cuts in Eq. (|10p . 

In this analysis, selection cuts are mass independent. As is seen, the total selection efficiency increases due to harder 
kinematic distributions with increasing This effect is, however, compensated by decreasing cross section of 

signal events. Table ITVl shows selection efficiencies for background samples. 
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FIG. 4: Transverse momentum distribution of the muon from r lepton decay in signal events in the Tfj,E™ %3S final state. 

Having obtained the number of signal and background events after all selection cuts, the signal significance is now 
calculated as 



Signal Significance 



N s 



(15) 



where Ns {Nb) is the final number of signal (total background) events after all selection cuts. Table fVl shows the 
signal significance for different miA,±\ values at 30 fb^ 1 . 



B. The n n ET sa Final State 



To search for this final state, the r jet identification is removed from the analysis. Instead, the focus is on muon 
identification and E™ lss distributions. Figure [§] shows the transverse momentum distributions of the muons. The 
muons in the events are required to satisfy the following requirements: 



p£ > 50 GeV, M < 2.5 



(16) 



An event is required to have exactly two muons passing the above requirements. As seen from Fig. [9j the soft muons 
in the Z+jets events which are due to the off-shell production of Z° bosons are removed by the above requirement 
and do not contribute to the signal selection. To study the number of jets in signal and background events, a jet 
reconstruction is performed. Figure [TOl shows the number of reconstructed jets which pass the following requirements: 



> 30 GeV, \r]\ < 2.5 



-"T , • (17) 

A jet is required to be separated enough from muons with px > 30 GeV in the event with the following requirement: 



ARn 



(jet,//.) 



> 0.4 



(18) 



10 




FIG. 5: Transverse energy distribution of reconstructed jets in signal and background events in the r/x_BJ* lss final state. All 
jets satisfying Eqs. [TT] and [T^] fill the histogram. Softer jets with Et < 30 GeV are discarded as the jet multiplicity strongly 
depends on the jet transverse energy and rapidly increases when soft jets are included, leading to long analysis time per event. 
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FIG. 6: Number of reconstructed jets in signal and background events in the r/ii5™ lss final state. A jet is considered and 
counted if it satisfies Eqs. llll and [T2l 
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FIG. 7: Number of tracks in the signal cone around the leading track of the t jet. 




FIG. 8: The distribution of the azimuthal angle between the muon and the r jet candidate. 
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80 GeV 


90 GeV 


110 GeV 


130 GeV 


Total cross section [fb] 


783 


521 


265 


150 


Number of events at 30 fb^ 1 


11745 


7815 


3975 


2250 


N Muons = 1 


5018(42.7%) 


3810(48.7%) 


2251(56.6%) 


1419(63%) 


N Jets = 1 


2029(40.4%) 


1585(41.6%) 


976(43.4%) 


647(45.6%) 


Isolation 


1375(67.7%) 


1101(69.4%) 


686(70.2%) 


468(72.4%) 


1- or 3-prong decay 


1306(95%) 


1045(94.9%) 


653(95.3%) 


447(95.5%) 


opposite charge 


1304(99.8%) 


1043(99.8%) 


653(99.9%) 


447(99.9%) 


total efficiency 


11.1% 


13.3% 


16.4% 


19.8% 


Expected events at 30/fe -1 


1304 


1043 


653 


447 



TABLE III: Selection efficiencies and remaining number of signal events after each cut in the TfiE™ 133 final state. Numbers in 
parentheses are relative efficiencies in percent with respect to the previous cut. Branching ratios have been taken into account 
in transition from the second to third row, thus the number of events at 30/fc -1 has been calculated as TV = a tot. x BR x Lumi., 
where a to t. is the value quoted in the second row, BR is the product of all branching ratios of decays taking into account 
the permutation factors, e.g. in this analysis it is 2 x BR((j> + — > tE™ ss ) X BR((j>~ — > /xi?™ ss ), and Lumi. is the integrated 
luminosity which is 30/6" 1 in this analysis. 



Process 


W+W~ 


ti 


W+jets 


Total cross section [pb] 


115.5 


878.7 


187.1 xlO 3 


Number of events at 30 fb~ 


577577 


4394086 


5.9 xlO 8 


N Muons = 1 


142084(24.6%) 


1.8c+06(40.9%) 


2.5e+07(4.2%) 


N Jets = 1 


60812(42.8%) 


81512(4.5%) 


l.le+07(44.5%) 


Isolation 


8483(13.9%) 


5533(6.8%) 


953682(8.6%) 


1- or 3-prong decay 


4807(56.7%) 


2545(46%) 


244810(25.7%) 


Opposite charge 


4618(96.1%) 


2241(88.1%) 


209655(85.6%) 


total efficiency 


0.8% 


0.05% 


0.03% 


Expected events at 30/6 -1 


4618 


2241 


209655 



TABLE IV: Selection efficiencies and remaining number of background events after each cut in the t^iE™ 38 final state. Numbers 
in parentheses are relative efficiencies in percent with respect to the previous cut. Branching ratios have been taken into account 
in transition from the second to third row. 

Thus, to suppress the background events with jets in the final state, it is required that no jet satisfying the above 
requirement exists in the event. Since Z+jets events could be a potential background for fifi final state, the invariant 
mass of the two muons are plotted as shown in Fig. QT] and a cut on that is applied as the following: 

Inv. Mass(/!i,/i 2 ) > 120 GeV . (19) 

In the next step the azimuthal angle between two muons is plotted as shown in Fig. 1121 In this analysis, we do not 
also apply any cut on A(f>^ l fi2 y As seen from Fig. 1121 the Z+jets events tend to show stronger back-to-back muon 
pairs. This effect can be due to the fact that muons in the Z+jets events come from a single mother particle. Further 
suppression of this background is possible by applying a cut on this distribution e.g. A(j)^ l fl2 - ) < 3. However, since 
the final signal significance exceeds 5<r, a study of the power of this cut and further optimization of this search is left 
for a full analysis including detector effects. Since the signal events are produced with opposite charges of muons, the 
opposite charge requirement is also applied to search for this final state: 

charge (Mi) + charge (Ma) = . (20) 





80 GeV 


90 GeV 


110 GeV 


130 GeV 


Signal significance 


2.8 


2.2 


1.4 


1 



TABLE V: Signal significance in the t/j,E™ ss final state for different m,^) hypotheses at 30 fb 
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Muon p T [GeV] 



FIG. 9: Muon transverse momentum distributions in signal and background events in the /x/i-E™ lss final state. This is an 
inclusive histogram i.e. any muon in the event satisfying kinematic cuts as in Eq. [TS]has filled the histogram. Histogram starts 
from 10 GeV on the X-axis to suppress the large peak of soft muons from off-shell Z+jets events. 
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FIG. 10: Number of reconstructed jets in signal and background events in the ^i/ii5™ ss final state. A jet is considered and 
counted if it passes kinematic requirements described in the text. 
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FIG. 11: Distribution of the invariant mass of the di-muon system in the /j,fj,E™ S3 final state. 
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FIG. 12: The distribution of the azimuthal angle between the two muons. 
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FIG. 13: Missing transverse energy distribution in the signal and background events in the fifiE™ t3S final state. 



Finally the E™ lss distribution is obtained using Eq. (J^TJ) . 

E™ ss = \J2p£\- (21) 

i 

The index i runs over all stable particles in the event within the volume defined by \r]\ < 3.5. This volume covers 
barrel and endcap calorimeters. The resulting distributions are shown in Fig. [13] and a cut on E™' lss is applied as the 
following: 

EJp zss > 50 GeV . (22) 

Table IVll shows the selection efficiencies and the remaining number of signal events after applying each cut for different 
m(0±) hypotheses. 

As can be seen the total selection efficiency increases with increasing m/^±y Table Ivnl shows selection efficiencies 
for background samples. In this final state, VF+jets background is more suppressed due to requiring two muons in 
the event. While one of the muons can be genuine from the W decay, the other has to come from heavy meson decays 
appearing as a muon inside a jet cone. These muons are expected to be better suppressed with more sophisticated 
reconstruction and isolation algorithms used in the LHC full simulation analyses. The Z+jets background is also 
expected to be under control with the cut on E™ lss and also the hard muon pair requirement. As can be seen from 
Fig. [9j the majority of these events tend to have either very soft muons or muons located around Et — 45 GeV 
(making the Z boson invariant mass peak). Thus, only a negligible fraction of these events appear with both muons 
harder than Et = 50 GeV. With the number of signal and background events as obtained in Tabs. IVII and IVII1 the 
signal significance is calculated as shown in Table IVIII1 As a rule of thumb, one can estimate the signal significance 
in the ee-E™ iss final state. For this search, the point B in Table U can be used. According to the values of g\ a which 
are equal for a = e, \i and t, the branching ratios of the (fy^ decays to different leptonic final states would be equal, 
i.e., Bi?((/) ± — > e ± Ni) — i?i?(0 ± — > H^Ni) — BR^*- -> t^Ni) ~ 1/3. Since the background samples have the same 
cross section times branching ratios for both final states, one can use the signal significance in mxE™ lss and scale the 
numerator of Eq. [15] by a factor of 4/9 which leads to a significance of 4.1cr for eeE™ S8 final state for m^±) = 80 
GeV at 30 /6 _1 . This estimation is based on the assumption that the kinematics of events in /j,/j,E t iiss and eei?™ ss 
final states are very similar. From another point of view this is a simplistic estimation as there is a difference between 
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80 GeV 


90 GeV 


110 GeV 


130 GeV 


Total cross section [fb] 


783 


521 


265 


150 


Number of events at 30 


5872 


3907 


1987 


1125 


N Muons = 2 


1348(23%) 


1080(28%) 


720(36%) 


494(44%) 


N Jets = 


961(71%) 


763(71%) 


494(69%) 


331(67%) 


Inv Mass(/i,^)>120 GeV 


895(93%) 


711(93%) 


463(94%) 


308(93%) 


Opposite charge 


895(100%) 


711(100%) 


463(100%) 


308(100%) 


rpmiss 


431(48%) 


393(55%) 


310(67%) 


231(75%) 


total efficiency 


7.3% 


10.1% 


15.6% 


20.6% 


Expected events at 30/6 - 


431 


393 


310 


231 



TABLE VI: Selection efficiencies and remaining number of signal events after each cut in the fifiE™ l3S final state. Numbers in 
parenthesis are relative efficiencies in percent with respect to the previous cut. Branching ratios have been taken into account 
in transition from the second to third row. 



Process 


W+W~ 


tt 


W+jets 


Z+jets 


Total cross section [pb] 


115.5 


878.7 


187.1 xlO 3 


258.9xl0 3 


Number of events at 30 fb" 1 


38713 


2.9xl0 5 


5.9xl0 8 


2.6xl0 8 


N Muons = 2 


4120(11%) 


52867(18%) 


2635(4xl0^ 4 %) 


752180(0.3%) 


N Jets = 


2925(71%) 


1705(3%) 


1219(46%) 


370599(49%) 


Inv Mass(^i, ^)>120 GeV 


2620(90%) 


1514(89%) 


1160(95%) 


220655(59%) 


Opposite charge 


2619(100%) 


1447(96%) 


1082(93%) 


220544(100%) 


Tpmiss 


730(28%) 


937(65%) 


295(27%) 


234(0.1%) 


total efficiency 


1.9% 


0.3% 


5xl0~ 5 % 


9xl0- 5 % 


Expected events at 30/6 -1 


730 


937 


295 


234 



TABLE VII: Selection efficiencies and remaining number of background events after each cut in the fifiE™ tas final state. 
Numbers in parentheses are relative efficiencies in percent with respect to the previous cut. Branching ratios have been taken 
into account in transition from the second to third row. 



muon and electron reconstruction and selection efficiencies in the detector. There can be also some under-estimation 
due to the fact that in the VF+jets events, the electron fake rate is higher than the muon fake rate. A detailed study 
of these features is beyond the scope of this paper and needs a full simulation of the detector. 



V. MEASURING THE MODEL PARAMETERS 



In this section, an approach for measuring model parameters is introduced and verified for the specific search 
described in this paper. 

Having observed N i, s at 30 fb" 1 , the expected number of signal events has to be extracted from the observed number 
of events by subtracting the contamination of the background and then dividing by the total selection efficiency of 
the signal events which is obtained from the simulation. In other words, 

N S = N ^-- Nb (23) 



m (0±) 


80 GeV 


90 GeV 


110 GeV 


130 GeV 


Signal significance 


9.2 


8.4 


6.6 


4.9 



TABLE VIII: Signal significance in /k/x_E™ !SS final state for different rn^±^ hypotheses. 



17 



Channel 


WW 


W+jets 


Z+jets 


tt 


Central Cross Section 


115.5 pb 


187.1 nb 


258.9 nb 


878.6 pb 


Ren. Scale x 2 


111.5(-3.5%) 


191.7(2.1%) 


263.1(1.6%) 


777.2(-11.5%) 


Ren. Scale -f- 2 


117.7(1.9%) 


181.8(-2.8%) 


252.2(-2.6%) 


977.3(11.2%) 


NNPDF2.0 


114.5(-0.9%) 


185.9(-0.6%) 


285.5(10.3%) 


917.6(4.4%) 


MSTW2008nnlo 


119.5(3.5%) 


198.1(5.9%) 


270.9(4.6%) 


861.9(-1.9%) 


CTEQ6.6 


113.3(-1.9%) 


196.3(4.9%) 


302.9(17%) 


841.7(-4.2%) 


a s =0.1127 


114.3(-1%) 


186.2(-0.5%) 


262.2(1.3%) 


813.6(-7.4%) 


a s =0.1207 


116.9(1.2%) 


187.9(0.4%) 


255.7(-1.2%) 


940.7(7%) 


Total uncertainty 


4.6% 


6.4% 


17% 


14% 



TABLE IX: Theoretical uncertainty of background processes excluding the uncertainty arising from the PDF set parameters 
errors. The central value of cross sections has been calculated using MRST2004 NNLO PDF set with a central value of 
a s = 0.1167. 



where es is the total selection efficiency of the signal, N$ is the expected number of signal events and Nb is the 
contamination in the observed sample which can be derived either from real data or simulation. In the latter case, 
the theoretical values of the backgrounds cross sections are used together with the selection efficiencies calculated 
from the simulation. Measuring N$ in two different final states would allow possibility of comparing two numbers 
from the two final states and estimating branching ratio of ± decays which in turn leads to the knowledge of model 
parameters such as gi a . 

In order to estimate the uncertainty with which N$ would be measured, the relative error on N$ is calculated as the 
following: 

AN S Ae s AN obs . AN B 

N S ~ e s ® Nf- ® Nf- (24) 

where © means square root of quadratic sum of all terms. Recall that N$ is the original number of events produced 
at the LHC, while N s eL is the number of selected signal events at the end of selection cuts; i.e., N s eL = egiVg. The 
relative uncertainty of the efficiency, Aes/es, is a quadratic sum of the main uncertainties; i.e., the jet energy scale 
uncertainty, missing transverse energy uncertainty and lepton reconstruction uncertaint y. As the best guesses for 
achievable uncertainties in LHC detectors, an uncertainty of 3% for jet energy scale Ref. [lj| an d missing transverse 
energy and 2% for lepton reconstruction and selection are assumed at an integrated luminosity of 30/6 -1 . The total 
selection efficiency could therefore have an uncertainty of 5% or better which is smaller than the uncertainty of the 
background in the sample due to the large amount of remaining background events. Suppose Nb is measured using 
theoretical cross section values and selection efficiencies from simulation. The uncertainty of Nb relative to signal is 
then 

JV| rf - ~ a L e B >Nf- ■ {2b) 

The uncertainty of the theoretical background cross sections originates from the choice of PDF set, the PDF set 
parameters errors, the value of as used in the PDF set and the renormalization and factorization scale uncertainty. 



Some of these uncertainties have been estimated for the tt events in Ref. [18J and lead to a ~ 14% uncertainty 
for this background which should be added to the uncertainty due to the PDF fit parameters. We also study 
the theoretical background uncertainties independently using MCFM package (excluding the error from PDF fit 
parameters which needs a dedicated detailed analysis). Results are shown in Tab. IIX1 The renormalisation scale is 
set to the characteristic mass in the event, i.e., the top quark, W boson and Z boson masses for tt, W+jets, WW and 
.Z+jets respectively. The variation of renormalisation scale is then performed around the central value. The central 
value of as according to PDG Ref. [l9| is 0.1176±0.002. The central value of this parameter chosen by different PDF 
sets is different, however, the 0.002 uncertainty is interpreted as one-cr or 68% uncertainty and a 95% uncertainty 
(~ 2a) equal to 0.004 is therefore used for the as variation. This leads to the values of 0.1127 and 0.1207 using the 
central as value adopted by MRST2004 NNLO PDF set. The total error has been calculated as proposed in Eq. [26] 

rp , , / (Ren. Scale t) 2 + (Ren. Scale lj* ~ f ~ ' (a s t) 2 + {as I) 2 

Total error = \ h (max PDF set error)^ H (26) 
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where (Ren. Scale f) an d (Ren. Scale i) are the observed error on the central cross sections when the renor- 
malisation scale is respectively set to twice and half the central value. Similarly, (as t) an d (c*s 4) are errors 
on central cross sections when as is set to respectively 0.1207 and 0.1127. For example, for the WW case, 
(Ren. Scale f) = -3.5 %, (Ren. Scale 1) = 1.9 %, (a s t) = 1-2 % and (a s 4) = -1 %■ The middle term in Eq. 
(|26|) is the maximum uncertainty observed due to a specific PDF set. As the error due to the PDF fit parameters 
can be of the order of the error due to the choice of PDF or as variation (e.g. cf. Ref. |20|) an average total 
error of ~ 10% is expected for the main background processes. It should be noted that here the errors due to the 
top quark, W and Z boson mass measurements were neglected although they could also contribute sizably to the 
total error. All arguments above imply that an average 10% error on the main backgrounds is a reasonable assumption. 

Taking 3% uncertainty for the LHC luminosity [21( and 10% uncertainty on background cross sections, for TO(^,±) = 
80 GeV, this term is estimated to be 

^prM - 58% ■ ( 2? ) 

Here only the result on the \x[i final state has been presented as the r/i final state suffers from the large background 
contamination. 

The second term on the right side of Eq. ([24| refers to the ratio of statistical uncertainty of the observed number 
of events to the number of signal events and can be written as 



N sel. N sel. ■ \ > 

This term amounts to 12% for the fifj, final states. Thus, the dominant contribution to ANs/Ns is the last term of 
Eq. (|24|) . As far as the uncertainties in the selection of events and the collider luminosity and cross section values 
remain the same, the last term, which is proportional to Nb/ N s eL , does not change with luminosity. An improvement 
is therefore possible only by acquiring a better knowledge of the uncertainties involved in the analysis or increasing 
the center of mass energy to obtain a higher ratio of signal to background. In principle, by using the data from the 
LHC itself, it would be possible to reduce the uncertainties. 



VI. AN ALTERNATIVE CHANNEL: THE ^fa + <j> ± 8 2 PRODUCTION 

In this section, we study the alternative channel, pp — > fa^fa + </ >± ^2- Remember that the gauge coupling of the 
neutral CP-odd and CP-even components of $ (fa and fa) to W^fa 1 are similar. On the other hand, for small values 
of the mixing (a ~ 0.01), $2 — fa- From the production point of view, 5% therefore behaves like fa so it is enough 
to study the first process (i.e., pp — > fa 1 fa) and double the number of events to account for the 4^82 production as 
well. In this analysis, we take m^ 2 = mj 2 = 90 GeV in accord with Ref. @ which imposes a lower limit of 90 GeV 
on the neutral Higgs boson mass. With the above assumptions, processes such as fa,82 — > 4 ± W T , fa — > 8 2 Z° or 
62 — > (f>2Z° will be forbidden and 4)2(82) will have only invisible decay modes. Since </>2 appears as missing energy in 
the detector, the final state to analyze is [iE™ %ss or tE™ %ss depending on the decay of ± to a muon or r. As before, 
we take BR(0 ± jt^iVi) = BR(0 ± r ± A^i) = 0.5. 

The simulation of these events is done by considering that fa, a CP-odd neutral scalar, would behave like A in 
MSSM and its production process would therefore resemble the H^A production in the framework of the MSSM. 
Within the MSSM, the CP-odd and CP-even Higgs particles are produced similarly. Since the #±,4° process is 
not produced by PYTHIA, we use H^H to produce the signal events. The coupling of to W^H is given by 
(e 2 /sin 2 6w) sin(/3 — a) so setting cos(/3 — a) = 0, the coupling will be the same as the coupling of ^ ± to W^&2 within 



our model. Moreover, in this analysis we set tan(/3) = 100 and obtain cross sections in agreement with Ref. [22j. To 
produce the same type of events, in production of H ± H° events, H° is forced to decay only to a muon pair and those 
muons are ignored in the event to have the same lepton multiplicity as in (fy^fa production. The missing transverse 
energy is also calculated by ignoring the H° products (muons) in the event so treating it as a source of invisible energy. 
Figure Q3] compares E™ lss as calculated only from the Ni in fa^ decay and the total E?p' lss including invisible energy 
from fa. Figure [15] shows the calculated total cross section of the signal for different ^ masses. This is the inclusive 
cross section although there is a little difference between the production rate of <p + and fa~ at the LHC. Since the 
final state is either ^E™ tss or tE™ iss , the main background for this signal would be single W boson production. All 
other backgrounds are expected to be suppressed by requirements on lepton multiplicity and central jet veto. Figures 
1161 1 1 71 1181 and [T5] compare the signal and background distributions of the muon transverse momentum, jet transverse 
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FIG. 14: Distribution of total E^" lss and its comparison with distribution of E™ lsa from <f> — > fiis. 
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FIG. 15: Total signal cross section in the ^ ± 02 channel as a function of ± mass. 
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FIG. 16: Muon transverse momentum distribution in <j> <j>2 events in fiE™ 11 " final state. 



energy, jet multiplicity and EJp iss respectively in the nE™ lss final state. Based on these distributions and since the 
background cross section is high and needs a large suppression factor, the following kinematic cuts are applied: 



Muon p T > 100 GeV, \rj\ < 2.5 



(29) 



E% lss > 100 GeV 



(30) 



Jet E T > 60 GeV, \r/\ < 2.5 



(31) 



the tE™ iss 



There should be no jet in the event with \xE™ lss final state passing the above requirement. Similarly Figs. [20l and 
compare the signal and background distributions of the jet transverse energy and jet multiplicity respectively in 
final state. The missing transverse energy distribution is not studied based on the same argument as 
mentioned in the T/j,EJp lss final state analysis of <p + 4>~ process. For this final state the same cut on the jet Et is 
applied. Since there is a r lepton in the event, the same r-jet identification as in the analysis of T[iE™ lss is applied. 
Similarly to the t/j,E™ iss case, the polarization dependent cuts are not applied. Based on selection requirements of 
the two final states, Tabs. |X] and IXII show selection efficiencies of signal and background events and the remaining 
number of events. Thus for m/0±) = 80 GeV, a signal significance of 9a and 4.6c is respectively obtained for /ii?™ ss 
and tE™ iss final states at 30 fb^ 1 integrated luminosity. Choosing the other point in the parameter space (the point 
B in Table|T|, one obtains Bi?((/) ± — > e ± Ni) = 1/3 so the signal cross section in the eE™ lss final state would be two 
third of that in the /j,EJp lss final state using point A in Table Q] The eE™ tss significance is therefore expected to be 
2/3x9 = 6er. Since the cross section of the signal decreases with increasing m^ij, higher cj^ masses would lead to 
lower signal significances. 



VII. DISCOVERY POTENTIAL AT THE 7 TEV RUN 



In this section, an estimation of signal significance at different final states is made by a simple rescaling of the signal 
and background cross sections to their corresponding values at the 7 TeV center of mass energy. Although a detailed 
study of kinematic distributions at this energy is needed to reach a concrete conclusion, a simple cross section rescaling 
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Jet E T [GeV] 



FIG. 17: Jet transverse energy distribution in </) ± 02 events in fiE™' sa final state. Jets satisfying Eq. [31] fill the histogram. 
Softer jets with Et < 60 GeV are skipped due to the same reason explained in caption of Fig. [S] 
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FIG. 18: Number of reconstructed jets for the <^) ± 02 events in iiE™ lss final state. A jet is considered and counted if it passes 
the kinematic requirements described in the text. 
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FIG. 19: Missing transverse energy distribution for the (j> y 2 even ts m the \iE™ %aa final state. 




FIG. 20: Jet transverse energy distribution for the 4> ±( t>2 events in the tE™ %S3 final state. Here also only jets satisfying Eq. [31] 
fill the histogram. Softer jets with Et < 60 GeV are skipped due to the same reason explained in caption of Fig. [S] 
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FIG. 21: Number of reconstructed jets for the <j> <\>i events in the tE™ laa final state. A jet is considered and counted if it 
passes the kinematic requirements described in the text. 



Final state 




Process 


Signal 


W+jets 


Total cross section [pb] 


1604 


187.3 xlO 3 


Number of events at 30 /6 _1 


24057 


5.9x10 s 


N Muons = 1 


3928(16%) 


1952900(0.3%) 


N Jets = 


3182(81%) 


356404(18%) 


E m,iss > 100 Ge y 


2652(83%) 


87319(24%) 


total efficiency 


11% 


0.015% 


Expected events at 


2652 


87319 



TABLE X: Selection efficiencies and remaining number of background events after each cut in the /i£"™ lss final state. Numbers 
in parentheses are relative efficiencies in percent with respect to the previous cut. Branching ratios have been taken into account 
in transition from the third to forth row. To calculate the signal, both (jy^fa and 4^82 have been taken into account. 



would give first feeling on how signal significance would decrease at lower center of mass energies. We choose the 7 
TeV center of mass energy to discuss the possibility of establishing the model introduced in this paper at the current 
LHC run. Table IXIII shows ratio of 7 TeV to 14 TeV cross section of signal and background processes. These ratios 
lead to the signal significance as shown in Table IXTlTl for different final states, provided that the data corresponding to 
30 fb~ x integrated luminosity is collected at the 7 TeV center of mass energy. This assumption might be in contrast 
to the current plan of the LHC machine, as there may be a switch to higher center of mass energies (10 TeV and 14 
TeV later on) before this amount of data is collected. If this phase finishes after collecting only 1 the expected 

signal significance will be l/y30 ~ 0.18 times smaller so there will be no discovery chance. However with 30 
as seen from Table IXTLlI searching for the <j) + (f>~ process, a 5 a discovery for the ii[iE™ iss final state for ± masses 
below ~ 100 GeV will be possible. A discovery for \xE™ tss final state of ± 02 + 4^ ^2 is also possible for mu±-\ = 80 
GeV. The needed integrated luminosity for this conclusion is 30 fb^ 1 . 
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Final state 


remiss 


Process 


Signal 


W+jets 


Total cross section [pb] 


1604 


187.3 xlO 3 


Number of events at 30 fb^ 1 


24057 


5.9x10 s 


N Jets = 1 


5117(21.3%) 


1.5e+07(2.6%) 


Isolation 


3058(59.8%) 


1.2e+06(7.7%) 


1- or 3-prong decay 


2900(94.8%) 


389380(33%) 


total efficiency 


12.05% 


0.066% 


Expected events at 30/6 -1 


2900 


389380 



TABLE XI: Selection efficiencies and remaining number of background events after each cut in the r_B™ lss final state. Numbers 
in parentheses are relative efficiencies in percent with respect to the previous cut. Branching ratios have been taken into account 
in transition from the third to forth row. To calculate the signal, both ± 02 and (^82 have been taken into account. 



Signal 


Channel 


Mass Point 


7 TeV to 14 TeV Ratio 




m w±) = 80 GeV 


0.4 


0+0" 


m w±) = 90 GeV 


0.39 




m w±) = HO GeV 


0.37 




m (<#.±) = 130 GeV 


0.35 


4^4)2 


m w±) = 80 GeV 


0.4 


Background 


Channel 


7 TeV to 14 TeV Ratio 


W + W~ 


0.38 


ti 


0.19 


W+jets 


0.49 


Z+jets 


0.6 



TABLE XII: Ratio of 7 TeV to 14 TeV cross sections of signal and background events. 

VIII. CONCLUSIONS 

In this paper, the possibility of discovering the SLIM model proposed in Q has been discussed. The focus has 
been on search for a signal at the LHC in two different categories based on the + _ and ± 02 + production 
processes. In the first category, two final states of r \i and \i /1 plus missing transverse energy have been analyzed in 
detail. The main background processes have been identified and a cut based search has been performed for the two 
final states separately. The search focused on 30/5" 1 integrated luminosity at LHC. It has been shown that the r /i 
and /x /i final states would have a signal significance of 2.8<r and 9.1a respectively for mu±) = 80 GeV. The ee final 
state has been estimated to have a significance of 4.1ct based on a simple rescaling of the fifi final state results. The 
above search allows to have a signal significance exceeding 5a in the /z/i final state for the ± masses below 130 GeV 
at 30/& -1 . In the second category, two final states of tE™ %ss and fiE™ lss were analyzed in detail. Results indicate 
that a significance of 4.6a and 9.8a can be obtained for the above final states respectively. The eE™ lss final state is 
expected to have a significance of 6.5a. A simple rescaling of the cross sections to their values at 7 TeV run of LHC 
shows that a discovery chance exists for miE™ lss final state of 4> + 4>~ and also for [i,E™ lss final state of 4> + (j)2 even 
at this center of mass energy, provided that data corresponding to 30 fb^ 1 integrated luminosity is collected before 
switching to higher machine energies. While the signal is observable at LHC, because of the uncertainties in the 
background estimation, measuring the new couplings of the model is more challenging. To make such a measurement 
a possibility, either these uncertainties have to be reduced or a future collider with a higher center of mass energy has 
to be employed. 

Apart from the signals that we have discussed in the paper (cj) + (f)~, ± 02 and (jr^S?), the model can have other 
observable effects such as the 52<t>2 production or invisible Higgs decay H — > SiSi 3]. Studying the prospect of 
detecting all these modes is beyond the scope of the present paper and will be presented elsewhere. 
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Channel 


Mass Point 


Signal significance 




m w±) = 80 GeV 


1.6 


4> + 4>- -> T/iEx iaa 


m w±) = 90 GeV 


1.2 




m w±) = 110 GeV 


0.7 




m {4>±) = 130 GeV 


0.5 




m w ±) = 80 GeV 


6.4 


(j> + <t>- -> fifiE^ 83 


m (tj>±) = 90 GeV 


5.7 




m w±) = 110 GeV 


4.2 




m w±) = 130 GeV 


3.0 




m (4l±) = 80 GeV 


2.6 


0*02 ->■ fiE™ iss 


m w±) = 80 GeV 


5.0 



TABLE XIII: Signal significance in different final states for the 7 TeV run, provided that 30 fb 1 data is collected at this 
energy before any switch to higher machine energies. 



Appendix A 

In this appendix, we derive constraints on the coupling <^ a from the neutrino mass matrix for the most economic 
case with only two right-handed neutrinos. As discussed in the text, in this case neutrino mass scheme is hierarchical. 

• Normal Hierarchical scheme: 

(m u ) a ^ = Upmns ■ Diag[0, y/ Am 2 sun , \J Am 2 atm e^} ■ U 



T 

PMNS 



where Upmns is the neutrino mixing matrix whose elements are known up to the subleading #13 effects. However, 
we do not know the value of the Majorana phase £. 

Using Eq. ©, we find that 



3 1 

where O is an arbitrary orthogonal matrix: T ■ O = 1 so 

O = 



(Al) 



cost/ suit 
— sin 9 cos ( 



(A2) 



In this case, we do not have any prediction for the coupling ratios independent of the arbitrary 9. 
• Inverted Hierarchical scheme: 



T 

PMNS ■ 



K) a p = Upmns ■ Diag[^/ Am 2 atm , \J Am 2 atm + Am 2 un e^,(j] ■ U 

From Eq. ([7)1. we therefore find 

9ia = £ j"(° T ' Diag[(Arr4j 1/4 , e' ?/2 (Am^ m + Am^) 1/4 ])« {UpmnsU , 
3 1 

where again O is an arbitrary orthogonal matrix as in Eq. IA2I It is noteworthy that regardless of the values 9 
and £, 



\9lr\ 2 _. \92r\ 2 



l + O{9 13l 23 -ir/4). 



|5imI 2 l32p| 2 

However the values of IgielVlgirl 2 or \g2e\ 2 /\g2r\ 2 depend on the unknown 9 and £. 



(A3) 
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